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P.E. Stallman and K. Itagaki 


INTRODUCTION 


It is the intent of this report to extend and improve the tables 
already available (CRREL Special Report 108, G.W. Aitken, 1966). The 
computational method is described with discussion of error, improved 
methods, and limitations. The tables are presented in two sections: 
the first for temperatures in the range -184°C to 0°c, the second for 
temperatures in the range 0°c to 100°C. The corresponding Fahrenheit 
temperatures are also included. 
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PROCEDURE 


The relationship of the emf produced by copper-constantan and the 
temperature is not well represented by a single function. Therefore, 
: for the purpose of generating these tables, two cubic equations with 
different sets of contants (A), Bj, C,) and (Adj, 32, Cz) were used. 
The cubic equation is 


Whe Cae RBAY 


es A,t +B,t? + ct? (1) 


where 7 = 1, 2, e = emf in millivolts, and t = temperature in Gy 


The constants A., 2. and C. in eq 1 were obtained by using the 
least squares polynotiial’curve Fit subroutine in the Dartmouth College 
GE 625 computer.* The input for the fitting function to obtain the 
constants Aj, 3,, Cy and A2,,B2, Cop was National Bureau of Standards 
Data (1955) from -190°C to 0°C and 0°C to 100°C respectively. All 
points given in the NBS tables in these ranges were used for the fits. 
The constants obtained from these fits were 


A, = 3.84971 x 1072 Ag = 3.84256 x 1072 
By = 4.62204 x 1075 By = 4.72065 x 1075 
Cy = -3.88666 x 1078 C2 = -3.76875 x 1078 


* This subroutine was written by L.C. Semprebon, Radiophysics Laboratory, 
Dartmouth College. 
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These constants were then incorporated into a computer program 
using the cosine method for the solution of the cubic equation (Weast 
et al., 1965). Thus, eq 1 is rearranged such that 


B A e 
$4242 4 4 22 <9: 
t? + C te + G t Cc 0 (2) : 
B =A aoe 
Define pP™aG? c’ r c’ 
let t=2-F (3) 
then zi+an+b=0 ; (4) © 
where a = + (3q - p*) and b = HF (2p3 - 9pq + 27r). at 
For the cosine solution @ is calculated by 
-(2p3 - 
cos ¢ = & SRE i ie ; (5) 
54 =a? a3 
27 27 


Then the particular solution for xz is calculated using 


2-2\2 (cos $+ 120°) . (6) 


Subs‘ ituting in eq 3, the temperature in °C is obtained. 


DISCUSSION 


As mentioned previously, a single function does not represent 
the emf-temperature relationship very well, and this was the reason 
for the two different fits. Even by using two fitting equations, there 
are some points that lie beyond a 0.5-microvolt confidence interval as 
shown in Figure 1. The main reason for this deviation from the NBS 
data is the rounding of the NBS values as well as the rounding of these 
tables to the nearest hundredth of a degree. Further information on 
rounding and improved methods using second degree Lagrange interpolation 
is given in the paper by Benedict and Ashby (1962). 


The tables given here are limited by the rounding discussed above 
and the physical limits of error produced by copper-constantan thermo- 
couples.* Considering these limitations, copper-constantan is preferred 
for work in the area covered by this table, in particular those tem-~ 
peratures below -50°C (Powell et al., 1962). 


*® ASTM Standards lists the limits of gtror for copper-constantan 
thermocouples between -75°F and 200°F as + 1-1/2 F for standard wire 
and + 3/4°F for special wire. For the range -300°F to -75°F the limits 
are given for special wire only as + 1% of the temperature measure. 
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curves from Nattonal Bureau of 
Standards data. 
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